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Abstract — We investigate the relationship between electronics 
and geometry in a nonplanar SiGe/Si resonant quantum cavity 
(RQC) subject to a magnetic field. The transfer matrix technique 
originally due to Usuki and coworkers [Phys. Rev. B, vol. 52, pp. 
8244-8255, 1995] has been modified to account for the nonzero 
local curvature of the RQC. Our results demonstrate that low- 
temperature ballistic magnetoconductance in nonplanar RQCs is 
highly sensitive to the changes in curvature for a wide range of 
magnetic field strengths. 

Index Terms — Si/SiGe nanomembranes, magnetotransport, 
curved nanostructures, geometric potential 



I. Introduction 

RECENT interest in nonplanar two-dimensional (NP2D) 
electron systems [1] has been driven by the emergence of 
reliable techniques for fabrication of flexible heterostructures, 
based on strained Si/SiGe [2], [3], [4] or GaAs [5], [6]. 
Curvature has become another degree of freedom available for 
manipulating electronic systems, which leads to novel basic 
physics and suggests applications of NP2D systems in na- 
noelectromechanical systems (NEMS) as ultrasensitive scales 
and sensors. The ability to design quantum devices using 
flexible substrates increases with our understanding of the role 
played by geometry in determining the electrical properties of 
these systems [7], [8], [9], [10]. Therefore, there is a need 
for versatile modeling tools capable of capturing the physics 
in nonplanar structures in a variety of geometries, biasing 
regimes, temperatures, and magnetic fields. To that end, we 
have made modifications to a transfer-matrix technique [11], 
[12], widely and successfully implemented at low magnetic 
fields and in planar 2D systems [13], [14], [15], in order to 
expand its functionality to include magnetic fields in non- 
planar geometries. 

In this paper, we present a calculation of low-temperature 
ballistic magnetoconductance in a nonplanar resonant quantum 
cavity (RQC) under low bias, based on the modified transfer- 
matrix method. The simulated 300x300 nm^ RQC (Fig. [B 
is formed on a (100) SiGe/Si hetero structure upon selective 
etching of the sacrificial oxide layer under the cavity's wings. 
Narrow (50 nm) quantum-point contacts connect the cavity to 
the wide leads (not depicted) at the source and drain ends; the 
rigid leads and the cavity spine are flat and tethered to the 
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Fig. 1. A nonplanar resonant quantum cavity on a SiGe/Si heterofilm, formed 
by selective etching of the sacrificial layer under the cavity's wings. Strain 
from the lattice mismatch provides the force that causes the wings to bow. 
Curvature is measured terms of a = W/{7tR), where W is the width of the 
cavity and R is the radius of curvature. 



substrate. The two-dimensional electron system is formed at 
the interface between the undoped strained Si and the relaxed, 
phosphorous-doped SiGe underneath (for example, epitaxially 
grown on SOI [4]). The structure can be topped with another 
layer of n-type doped SiGe to aid in the electronic confinement 
[4] ; such SiGe/Si modulation doped heterostructures have been 
shown to produce good electronic confinement both before and 
after selective etching of the sacrificial oxide layer. But unlike 
the strain-compensated structure of Ref . [4] that remained flat 
upon underetching, our structure releases the wings to curl up 
into a partial cylinder, similar to what happens in the inverted 
Si/SiGe structure described in Ref. [3] that curled downward. 
Curvature is measured in terms of a parameter a = W/ (ttR), 
where W is the width of the cavity, and R is the radius of 
curvature. 

In Sec. ini we introduce the modified transfer-matrix 
method, applicable for low-bias transport in nonplanar struc- 
tures under magnetic field. Simulation results illustrating the 
effects of curvature on transport are presented in Sec. [nil and 
concluding remarks are given in Sec. [iVl 

II. The Modified Transfer-Matrix Method 

The transfer-matrix method due to Usuki et al [11] is 
an iterative, stable routine used to solve the discretized two- 
dimensional Schrodinger equation. The entire numerical do- 
main is divided into slices perpendicular to the current flow. A 
transfer matrix is constructed for each slice and, after iterating 
slice-by- slice in the direction of the current flow, the space- 
resolved electron density and the structure's transmission 
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coefficient are obtained. In the absence of a magnetic field, 
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(1) 



^n,n-i = diag[-£:o], 



where ^/^^ is the wavefunction of the n-th sHce, Tn is the 
transfer matrix from sHce n — 1 to sHce n, Hn,n±i are the 
hopping Hamiltonian matrices that connect nearest-neighbor 
sHces, and Hn is the finite-difference representation of the 
Hamiltonian of the n-th slice. All these Hamiltonians are 
matrices of dimension M, where M is the number of grid- 
points in the direction perpendicular to the current flow, while 
ipn is an M-column. Moreover, Eq = fi^/2m*a^ is the 
nearest-neighbor hopping energy, where a is the mesh spacing 
(assumed constant throughout the structure), and m* is the 
electron effective mass. Since Si is biaxially strained in a 
{100} plane, among its 6 equivalent valleys the two whose 
heavy mass is along the growth direction come down in energy 
with respect to the other four (the so-called valley splitting 
into two lower A2 and four higher A4 valleys [16]). The 
electrons in A2 valleys respond to in-plane bias with the 
light mass m* = 0.19mo (mo is the free-electron rest mass), 
which is the effective mass we use in the calculation. The 
subbands emerging from the A2 valleys upon confinement 
are considered twofold degenerate (although the magnetic field 
can introduce further splitting [16]). We will assume that only 
the lowest (twofold degenerate) A2 subband is occupied. The 
total electron sheet density in the leads, 2N, and the Fermi 
level Ef measured with respect to the subband bottom are 
therefore connected by 

tu'^Ef 



N : 



(2) 



In order to calculate conductance though the structure, the 
transmission coefficient for each propagating mode injected 
from the source is calculated. At the source end of the device, 
there is a wide lead (not depicted in Fig. 1), where many plane 
waves are injected with unit amplitude into the device's active 
region (cavity plus quantum-point contacts); only some make 
it all the way through to the drain lead, with their amplitude 
significantly diminished. Plane wave injection at the source 
end is achieved by solving Eq. ([T]) forn = and tpi = ipoe'^^^ 
(phase shift between two adjacent slices will be ka when a 
plane wave is injected, where k is the wave number), so 



Aka 



To 



^1 



(3) 



Since the elements of the transfer matrix To are real [Eq. ([T])] 
in the absence of magnetic field, the forward and backward 
propagating modes occur in pairs: of the 2M modes obtained, 
M are forward and M are backward modes; for each direction, 
a small number of those M modes are propagating and the 
rest are evanescent, and the latter do not contribute to the 
conductance. Only the forward propagating modes enter into 
the conductance calculations, and thus must be selected from 



the collection of e*^^ eigenvalues. This is a straightforward 
process in the absence of magnetic field: all propagating 
modes will have real k values, and can be identified as 
|gz/ca| _ ^ Forward modes will have positive k values, and 
thus imag(e*^^) > 0. Note that the grid spacing a must be 
sufficiently small to ensure that all ka < n. 



A. Modifications for magnetic fields 

In the presence of a magnetic field, a gauge-invariant 
form of the Schrodinger equation with the magnetic field 
incorporated can be obtained through the Peierls substitution 
[17]. This substitution has the effect of producing the (p+eA)^ 
term in the full Hamiltonian, which in itself is not suitable for 
implementation in transfer-matrix approaches. 

At zero magnetic field, one can obtain the matrix ele- 
ments Hpq of a tight-binding Hamiltonian by discretizing 
the Schrodinger equation. The effect of the magnetic field 
is introduced by modifying each Hamiltonian matrix element 
Hpq by an appropriate Peierls phase factor, containing the 
integral of the magnetic vector potential along a straight line 
connecting meshpoints p and q [18], [19], [20]: 



pq 



H. 



pq 



A-d£ 



(4) 



In the Landau gauge A — (—5^,0,0), the presence of the 
magnetic field in the transfer matrix method is felt only in 



f ^n,n+i) = -£^0 exp (z27rj<l>/<l>o) , j = 1, . . . , 



M, 



-'-'n,n+l> 



(5) 



where ^ = Bo? is the flux through a meshcell, and is the 
magnetic flux quantum. 

When the magnetic field is present, the matrix elements 
of To become complex, and the forward and backward prop- 
agating modes are no longer mutual complex conjugates. 
Since the essence of the transfer-matrix technique is proper 
initialization of forward propagating modes, the selection 
process becomes more complicated. Namely, the eigenvalues 
ex.p{ika) the phase shifts associated with the modes all get 
an extra positive phase shift proportional to ^/^o, which 
is roughly an average of all the Peierls' phases within one 
slice [12]. Consequently, some modes that are actually back- 
ward propagating may erroneously be identified as forward 
propagating (see the middle panel of Fig 121). With increasing 
magnetic field, one can always first try to decrease a, the 
meshsize, to minimize the phase shift, but one quickly runs 
into problems with computer memory (number of meshpoints 
in the structure and consequently the dimensions of matrices 
dealt with become too large). Instead, we have found that a 
simple shift of the coordinate origin - placing the origin at 
the center of the slice rather than the edge (Fig|2|)- keeps this 
extra phase shift virtually negligible up to very high magnetic 
fields. Certainly, this is not a fundamental physical issue, but 
rather a particularity of the transfer matrix technique. 
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Fig. 2. Absolute values and imaginary parts of all modes for a 300 nm 
wire at = 4 X 10^1 cm -2 and B = (top graph) and B = 0.4T (middle 
and bottom graphs). The Peierls phase skews the imaginary values of all 
eigenvalues and erroneously reports an extra forward propagating mode (red 
box) under this moderate magnetic field (middle graph). This skewing effect 
is eliminated by taking f A • d£ from the center of each stage (bottom graph) 
rather than the edge. 

B. Modifications for nonplanar geometries 

To investigate the effects of curvature combined with the 
presence of a magnetic field, the transfer matrix implementa- 
tion must be further modified. 

In the absence of a magnetic field, the motion of a particle 
confined to a curved surface can be described in the surface's 
natural coordinates, but with the addition of an extra attractive 
geometric potential [21]. The geometric potential Vg depends 
on the mean and Gaussian curvatures of the surface at a given 
point. For example, in the case of a cylindrical surface of 
radius R, the geometric potential is a constant and is given 
by = —h'^/8m*R'^. Indeed, the transfer matrix method as 
described in Seclllis perfectly applicable to other geometries, 
provided that one accounts for Vg, obeys the structure's 
modified metric when constructing the kinetic energy operator 
(details can be found in DaCosta's paper [21]), and discretizes 
the Schrodinger equation in the surface's natural coordinates. 

When a magnetic field is added to the motion over a 
nonplanar surface, the line integral J A ■ d£ that is used in 
modifying the Hamiltonian matrix elements Hpq must now be 
taken over a geodesic between p and q rather than along a 
straight line. This implies that one first needs to calculate the 
equations for the geodesies for the surface, and then integrate 
over those to obtain the Peierls phase. 

III. Results 

The aforementioned modifications were made to the transfer 
matrix method in order to observe the effects of curvature on 
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Fig. 3. Conductance variation with curvature at B = 0. At particular electron 
densities in the leads, the conductance is highly dependent on the curvature 
of the cavity even in the absence of a magnetic field. This phenomenon can 
only be attributed to the geometric potential. Insets: colorized plots of the 
magnitude squared of the wavefunction at the Fermi energy, for the specific 
values of a shown. 

the electronic properties of the nonplanar RQC depicted in 
Fig. 1. Our system is not fully cylindrical (it has cylindrical 
wings and a flat spine) and is therefore more complicated 
than a cylinder due to the varying geometric potential: Vg 
is constant and attractive in the wings, and zero in the 
spine. (An interested reader can find the analytical solution 
to magnetotransport on a cylinder with very small curvature 
in Ref. [5].) 

In Figs. 3-6, conductance variations in the structure (in units 
of Go = 2e^//i, the conductance quantum) are presented. All 
the results are obtained in the low-bias regime as well as at 
low temperatures, so the electron transport occurs at the Fermi 
level, specified by N, the electron sheet density (per A2 valley) 
in the source and drain leads [Eq. Q] . Conductance calculated 
is also per A 2 valley. Conductance plots are accompanied by 
color plots of the magnitude squared of the wavefunction at the 
Fermi energy, with low/moderate/high values corresponding 
to the usual blue/yellow/red color scheme. (This is the same 
wavefunction that determines the transmission coefficient. It 
involves many plane waves being injected with unit amplitude 
from the source lead into the structure's active region, as 
described in Section II. ) 

A. Zero Magnetic Field 

The geometric potential defined in III-BI is quite weak at the 
dimensions of our RQC, and its effect often goes undetected. 
At certain sheet electron densities TV, however, it alone is 
responsible for significant variation in conductance. Hence, 
conductance can be modulated by changing the curvature 
even if in the absence of a magnetic field, as shown in Fig. 
[51 At lower curvatures, most of the plane waves incident 
from the source lead (top) are transmitted to the drain [note 
the high probability density in the drain-end quantum-point 
contact (right-hand- side inset)]. As the curvature increases, the 
coupling between the source and drain decreases, leading to a 
decrease in conductance. 
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Fig. 4. Conductance variation with curvature at a low magnetic field (B = 
0.4T). Insets: colorized plots of the magnitude squared of the wavefunction at 
the Fermi energy, for the specific values of a shown. Patterns observed in the 
probability density are altered by changes in curvature, significantly affecting 
the conductance. 



B. Low Magnetic Field 

At low magnetic fields, the resonance patterns which form 
in the cavity on account of the magnetic field [14], [15] can 
be altered by cavity curvature. These patterns determine the 
transmission through the cavity (see Fig. |4]), and hence alter 
the conductance as well. 

C. High Magnetic Field 

At high magnetic fields, the results of our calculations con- 
firm what is known about the physics of high-field magneto- 
transport in planar structures, such as the formation edge states 
and depopulation of Landau levels. In addition, novel effects 
of curvature on current conduction in the strong quantum Hall 
regime are also observed (see FigQ. At high fields, the edge 
states (top panel of Fig. [5]) become wider and drift away from 
the edge of the cavity as the cavity rolls up (its curvature 
increases). As the cavity takes on a distinctly cylindrical shape, 
a mid-cavity state is formed as a result of opposing fluxes 
through the top and bottom portions of the wing. 

At certain magnetic fields (Fig. [6]), even minuscule changes 
in the curvature cause significant changes in the conductance. 
The conductance will remain sensitive to the changes in cur- 
vature in the high-field regime up to extremely high magnetic 
fields (B > ST in our particular example), after which the 
edge states become increasingly difficult to disturb by changes 
in curvature. 

IV. Conclusion 

We have developed a modified transfer-matrix method that 
enables low-bias, low-temperature transport calculations for 
nonplanar two-dimensional electron systems in the presence 
of magnetic fields. We have shown that the interplay be- 
tween geometry and electronic properties of flexible SiGe/Si 
heterostructures (such as the simulated cylindrical resonant 
quantum cavity) is substantial. In particular, our results demon- 
strate that conductance through a cylindrical RQC at a given 




Fig. 5. Magnitude squared of the wavefunction at the Fermi energy in the 
nonplanar RQC as a function of curvature at a high magnetic field (B = 
12. 5T). N = 4 X 10^^cm~2. At low curvature (a = 0.4, top panel) the 
edge state is clearly visible. As curvature increases (a = 0.5, second panel), 
the edge state migrates away from the edge of the cavity and broadens. As 
the cavity takes on a distinctly cylindrical shape (a=1.5, third panel and 1.7, 
bottom panel), it exhibits a mid-cavity state, created by opposing magnetic 
fluxi ility 
denf the 
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Fig. 6. Conductance variation with curvature at a high magnetic field (B = 
6T). Insets: colorized plots of the magnitude squared of the wavefunction at 
the Fermi energy, for the specific values of a shown. At certain magnetic 
fields, such as the one presented here, the cavity allows edge states which 
are prohibited in the narrow quantum-point contact. The shape of the cavity 
determines how these edge states couple to the modes in the contact; therefore, 
conductance again varies with curvature. 
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carrier density in the leads is highly sensitive to the changes in 
curvature for a wide range of magnetic field strengths, which 
suggests possible applications of nonplanar structures as nano- 
electromechanical systems (NEMS) and sensors (e.g., a dilute 
chemical presence that affects the shape of the film could be 
detected electronically), at least at very low temperatures. At 
higher temperatures, scattering due to phonons must not be 
disregarded, and energies within a few ksT around the Fermi 
level will contribute to the conductance, which may require 
including several A 2 subbands as well as some low-lying A4 
subbands. Consequently, smearing of the resonant features in 
G\s a (Figs. 3, 4, and 6) ought to be expected. Exploration of 
non-cryogenic transport in nonplanar RQCs will be the topic 
of a subsequent publication. 
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